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Abstract: We discuss the importance of event-by-event fluctuations in relativistic heavy-ion collisions. We 
present results of multiplicity fluctuations of photons and charged particles produced in 158 A GeV Pb-nPb 
collisions at CERN-SPS. Multiplicity fluctuations agree with participant model calculations indicating the 
absence of critical fluctuation. Localized fluctuations in charged particles to photons have been analyzed 
and compared to various mixed events. The results show localized fluctuations in both charged particles and 
photons. However, no correlated DCC-like fluctuation is observed. 

The primary goal of ultra-relativistic heavy-ion collisions is to investigate the properties of matter at high energy 
density and high temperature. These studies provide information on the dynamics of multiparticle production mech- 
anism and the phase transition from hadronic matter to Quark-Gluon plasma (QGP). One of the major advantages of 



heavy-ion collisions at high energies is that most of the observables can be studied in an event-by-event (E-by-E) basis 
because of the production of large number of particles in every event. This is particularly important as one expects 
large fluctuations in several observables associated with the phase transition. 

The subject of E-by-E fluctuations has recently gained considerable interest primarily motivated by the near perfect 
Gaussian distributions of mean transverse momentum and particle ratios [1] obtained by the NA49 experiment at 
CERN-SPS. As we will see later in this manuscript the multiphcity distributions for narrow centrahty bins are also near 
perfect Gaussians. The variance of the width of these distributions contain information about the reaction mechanism 
as well as the nuclear geometry. The task here is to distinguish between statistical fluctuations and those which have 
dynamical origin. Several methods have been put forward which suggest ways to infer about the presence of dynamical 
fluctuations which have new physics origin. These are done by evaluating the effect of nuclear geometry, influence of 
hadronic resonances and rescattering, and contributions from Bose-Einstein correlations [2, 3, 4, 5, 6]. Primary interest 
would be to identify critical fluctuations associated with the QGP phase transition. 

Another interesting phenomenon is the formation of disoriented chiral condensates (DCC) [7], which is a conse- 
quence of chiral phase transition. It gives rise to large E-by-E fluctuations in the number of charged particles and 
photons in localized phase space. The distributions of correlation functions using charged particles and photons have 
been shown to be near perfect Gaussians in the absence of DCC. The presence of DCC tends to make the distributions 
wider as has been shown in ref. [8]. 

In order to infer about the presence of non-statistical fluctuations, one has to compare the distributions obtained 
from experimental data to known models which incorporate all the known phenomena. An alternate or may be 
comphmentary procedure to probe these fluctuations in a model independent manner would be to compare data 
distributions to those of mixed events generated from the data. Once properly understood, mixed events provide the 
most unambiguous means to probe any new physics. In this manuscript we present results from Pb+Pb collisions at 
158 • AGeV taken at CERN-SPS by the WA98 experimental setup. The results will be compared to a participant model 
and also to several different types of mixed events. 

In the WA98 experiment, the emphasis has been on high precision, simultaneous detection of both photons and 
hadrons. Photon multiplicities (A'y-iike) were measured by the preshower photon multiplicity detector (PMD) placed 
at 21.4 meters from the target and covering an r| range of 2.9 —4.2. Charged particle hits (A'ch) were counted using 
a circular sihcon pad detector (SPMD) located 32.8 cm from the target with a coverage of 2.35 < rj < 3.75. The 
centrahty of the interaction is determined by the total transverse energy (Ef) measured in the mid rapidity calorimeter. 
The centralities are expressed as fractions of the measured cross section from the minimum bias Ej distributions. 

Figure 1 shows the minimum bias distributions of Wy-iike clusters and A'ch- Superimposed to these are the distri- 
butions corresponding to the centrahty cuts of top 1%, 2% and 5% of the minimum bias Ej distributions. These 
distributions turn out to be near perfect Gaussians with X^/ndf close to unity. The solid lines show the Gaussian fits 
to these curves. It has been observed that the distributions deviate from Gaussians if the centrahty bin is made broader 
than — 5%. In order to study the centrahty dependence of fluctuations we have chosen 2% bins in centrahty, viz-, 
0-2%, 2-4%,. ..,58-60%,.. etc. The resulting multiplicity distributions are near perfect Gaussians over large centrahty 
range, from peripheral to central collisions. 

If the distribution of a quantity, X, is Gaussian, then the amount of fluctuation may be defined as: 

cox = ^, (1) 

where Ox is the variance of the distribution, and {X) denotes the mean value of X. In order to compare the fluctuation 
for peripheral to central colhsions, we have calculated the relative fluctuation defined as (Ox per number of participant 
nuclei, A^pan (obtained from VENUS event generator). 

Relative multiplicity fluctuations of photons (cOy/A'part) and charged particles (C0ch/A'part)> as functions of A^part are 
shown in left and right panels of Figure 2, respectively. Since the measured data for photons contain contaminations 
from charged particles, these have been corrected to obtain distributions for photons only by using the efficiency and 
purity of the measured photon-like clusters [9]. The resulting data values have been compared with simple calculations 
based on participant model [6, 10] and also those obtained from VENUS event generator. A reasonable agreement have 
been seen with the participant model as well as VENUS indicating the absence of any critical fluctuations. 

Theoretical predictions suggest that the isospin fluctuations caused by DCC would produce clusters of coherent 
pions in phase space forming domains localized in phase space. We present results based on E-by-E fluctuation in the 
relative number of charged particles and photons. Details may be obtained from Ref. [11]. Two different methods are 
used for DCC search. The first method deals with the correlation of Atj. uke and A^ch and the second method is based on 
the discrete wavelet techniques (DWT). 
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FIGURE 1. Minimum bias distributions of (a) y-like clusters and (b) cliarged particle multiplicity produced in Pb induced 
reactions at 158-A GeV on Pb. Superimposed are the multiplicity distributions for the top 1%, 2%, and 5% most central events. 
These distributions turn out to be near perfect Gaussians as seen by the soUd fitted lines. 




FIGURE 2. The relative fluctuations of (a) photons and (b) the charged particles compared to calculations from a participant 
model and VENUS. The photon data have been corrected from AEy-iite to Ny using the efficiency and purity of the measured photon 
sample. Reasonable agreement has been seen between data and models. 



The correlation of A'y-iike and A'ch has been studied in smaller (|)-segments by dividing the (t)-space into 2,4,8 and 
16 bins. A common correlation axis (Z) has been obtained by fitting the distribution of the correlations with a second 
order polynomial. The distribution of the closest distance of the data points to the correlation axis represents the 
relative fluctuation of Ny-uke and A'ch for any given ^ bin. The final distributions (Sz) for 4 and 8 bins are shown in the 
left- top panel of Figure 3. Superimposed are the distributions from GEANT simulated events and a particular set of 
mixed events, called Ml. 

We have also used the DWT method where the analysis is performed by dividing the azimuthal space of PMD and 
SPMD into 2-i bins where j is called the scale of the division. The input to the DWT analysis is a spectrum of the 
function, Ny^/ {Nyi^ +Nch), at the highest resolution scale, 7max(= 5). The output consists of a set of wavelet or the 



TABLE 1. Four different types of mixed events are constructed by combinations of PMD and SPMD hits. 
The last column provides the usefulness of each of the mixed events to probe different type of fluctuation. 
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M3-ch 


Mix hits 

Unaltered hits 
Unaltered hits 
Mix hits 


Mix hits 
Unaltered hits 

Mix hits 
Unaltered hits 


Correlated + Individual 

Correlated 
Nj only 
Nch only 



so called father function coefficients (FFC) at each scale, from j = 1 to j = jmax — 1- The FFC distributions at 4 and 
8 bins in (|) are shown in the left-bottom panel of Figure 3. The Sz and FFC distributions are Gaussian in nature in the 
absence of DCC. The presence of events with DCC domains of a particular bin size in (|) would result in a broader 
distribution of Sz and FFC distributions compared to normal events [8, 1 1]. 

The presence of DCC-like fluctuations in the experimental data can be inferred by comparing the data distributions 
to different sets of mixed events. Properly constructed mixed events preserve all the detector effects while removing 
correlations. Four sets of mixed events are generated by maintaining the global (bin 1) A'yiike'^ch correlations as in 
the real data. Table 1 summarizes the construction of mixed events and the type of fluctuation these can probe. These 
behavior has been understood by simulating the effect of DCC and analyzing the simulated data in a manner identical 
to the real data. 
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FIGURE 3. The Sz and FFC distributions for 4 and 8 divisions in (|) are shown in the left panel. The experimental data. Ml and 
V-i-G events are shown by solid circles, soUd histograms and dashed histograms, respectively. The statistics for data and mixed 
events are the same, whereas the distribution for the V-i-G events is normaUzed to the number of data events. The right panel shows 
the root mean square (rms) deviations of the Sz and FFC distributions for various divisions in the azimuthal angle. 

The rms deviations of the Sz and FFC distributions as a function of the number of bins in azimuth are shown for 
experimental data, four types of mixed events and Geant simulated data in the right panel of Figure 3. The statistical 
errors are small and the systematic errors are shown in the figure. 



First we compare the rms deviations obtained from the data to those of Ml mixed events. The data points are much 
larger compared to those of Ml for both Sz and FFC distributions. For 2,4 and 8 bins the values of the rms deviations 



of data are 2.5a, 3.0a and 2.4a larger than those of the Ml events. Similarly for 4 and 8 bins the values of the FFC 
rms deviations of data are 3.7o and 2.8o larger than those of the Ml events. These indicate the presence of localized 
non-statistical fluctuations. The observed effect may arise because of (a) E-by-E correlated Afy_iike-A^ch fluctuations, 
(b) fluctuations in A'y-iike, and (c) fluctuations in Nch- The source of the fluctuation can be obtained by comparing the 
data with M2 and M3 types of events. 

Next we compare data rms deviations to those of M2. It is observed that in all cases the results from M2 match 
completely with data. This implies the absence of any E-by-E correlated (DCC-like) fluctuations and allows us only 
to put an upper limit on DCC production [11]. 

Presence of individual fluctuations are to be inferred by comparing data with M3 type of events. The M3 type of 
events are found to be similar to each other within quoted errors and he between Ml and M2. This indicates presence 
of localized independent fluctuations in both A^y-uke and A^ch- 

In summary, we have presented multiplicity distributions and multiplicity fluctuations of photons and charged 
particles produced in 158v4 GeV Pb+Pb collisions at the CERN-SPS. Multiplicity fluctuations from data agree 
reasonably well with participant model calculations over a wide range of centrahty, from peripheral to central 
collisions. This indicates the absence of any critical fluctuations indicative of phase transition at SPS energies. 
Localized fluctuations in photons and charged particles have been studied in order to search for formation of domains 
of DCC. The results from correlation and DWT analysis are compared to those of mixed events generated from data. 
The results provide model independent evidence for non-statistical fluctuations in the data for ^ interval between 
45 and 90 degrees and one unit in pseudorapidity. This is seen to be due to individual fluctuations in both photons 
and charged particles. No significant E-by-E correlated fluctuations are observed, contrary to expectations for a DCC 
effect. 
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